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Analytical models describing the motion of colloidal particles in given velocity 
fields are presented. In addition to local approaches, leading to well known master 
equations such as the Langevin and the Fokker-Planck equations, a global description 
based on path integration is reviewed. This shows that under very broad conditions, 
during its evolution a dissipative system tends to minimize its energy dissipation 
in such a way to keep constant the Hamiltonian time rate, equal to the difference 
between the flux-based and the force-based Rayleigh dissipation functions. At steady 
state, the Hamiltonian time rate is maximized, leading to a minimum resistance 
principle. In the unsteady case, we consider the relaxation to equilibrium of harmonic 
oscillators and the motion of a Brownian particle in shear flow, obtaining results that 
coincide with the solution of the Fokker-Planck and the Langevin equations. 

I. INTRODUCTION 

In classical thermodynamics, an isolated system tends to assume its stable equilibrium 
configuration which, according to the second law, is unique and corresponds to the state x in 
which the entropy functional S (x) is maximized. As in classical thermodynamics any system 
has an infinite number of degrees of freedom, the equilibrium state is infinitely more probable 
than any other state. On the other hand, finite systems fluctuate, as all configurations have 
a finite probability to occur, equal to the Gibbs distribution, C exp [S (x) /k \, where k is 
Boltzmann’s constant, and C is a normalization factor. 

Most of times, however, we deal with evolving dissipating systems, whose trajectories 
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are described through appropriate master equations; in particular, we are interested either 
in systems that are maintained at steady state through imposed boundary conditions, or 
systems that are relaxing towards their state of stable equilibrium. Now, the evolution of 
conservative systems can be described either through a mechanistic approach, by integrating 
Newton’s equations of motion, or through the variational principle of least action, leading 
to the Lagrangian and Hamiltonian formulations of classical mechanics. So, it is natural to 
look for a global description of irreversible processes as well, that is finding some universal 
function that could play the role of the action and the Lagrangian functionals in analytical 
mechanics. Naturally, just like the principle of least action and Lagrangian mechanics are 
equivalent to Newton’s equation of motion, this novel variational approach should constitute 
a way to describe irreversible processes that is totally equivalent, although alternative, to 
solving the usual master equations. 

The first example of a similar variational law is due to Kirchhoff-, who in the mid nine¬ 
teenth century demonstrated that the steady state transport of electric charge obeys the 
principle of least dissipation of energy. This principle was generalized one hundred year 
later by Prigogine^, stating that under given boundary conditions, a system will tend to 
a state, in the eventual stationary process, that has a minimum of entropy production. 
De Groot and Mazur- applied this principle to heat conduction, basically using the same 
approach as Kirchhoff’ in electrostatics, stressing however that the principle of minimum 
entropy production is true only when the corresponding thermodynamic flux (i.e. the heat 
flux, Jq ) is proportional to its conjugated force (i.e. the gradient of the inverse temper¬ 
ature, VT” 1 ) through a phenomenological coefficient that is constant, i.e., independent of 
temperature^. A discussion about the history of this principle can be found in Jaynes^ and 
in M filler and Weiss^. 

An apparently opposite principle was proposed by Ziegler 7 -, stating that when the thermo¬ 
dynamic forces, F, are assigned, then the actual fluxes, J, maximize the entropy production 
rate, S = J • F. This principle generalizes the maximum (plastic) dissipation theorem of 
Mises, Taylor and Hill 8 to all nonequilibrium thermodynamics. In particular, when thermo¬ 
dynamic fluxes and forces are linearly related to one another, i.e. J = L • F, with L ik = L kl 
denoting a constant generalized conductivity, then this theorem reduces to the well known 
principle, proposed by Onsager- in 1931, stating that when the thermodynamic forces, F, 
are assigned, the difference S — T j between the entropy production rate and the flux-based 
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Rayleigh’s energy dissipation function, \1 >j = | (J • L _1 • J), is maximized. Vice versa, when 
thermodynamic fluxes are assigned, the difference S — ^ p between the entropy production 
rate and the force-based Rayleigh’s energy dissipation function, fyp = | (F • L ■ F), is max¬ 
imized. Finally, for generic constraints, Onsager and Machlup^ and GyarmatiR- found that 
the difference S — \l/ j — \l/i? between the entropy production rate and the sum of the two the 
energy dissipation functions, is maximized. This extremum property is generally referred 
to as the maximum entropy production (MEP) principle; its implications are reviewed in 
detail by Martyushev and Selezne v 1213 and by Verhas^, showing, in particular, how the 
parabolic equations appearing in heat and mass transport can be derived from it. Finally, 
Berett a 1516 extended the MET principle to far-from-equilibrium conditions, reformulating it 
as the steepest entropy ascent (SEA) model. 

In this article, without trying to solve this fundamental problem in its entirety, we con¬ 
centrate on a particular case, namely the time evolution of a fluctuating system, subjected 
to conservative forces. Apart from its simplicity, the great advantage of this problem is that 
for linear forces it reduces to the Ornstein-Uhlcnbeck problem—, therefore providing a way 
to check the validity of our results. The novelty of the approach presented here is that we 
use the path integral formulation, thus providing an alternate point of view which might 
help to solve similar problems in the future. 


II. DESCRIPTION OF THE MODEL 

Consider a diffusing Brownian particle, subjected to a conservative force F = — V0O; 
satisfying the Langevin equation, 


C-(x-V) = f, (1) 

where the dot indicates time derivative, £ is the resistance dyadic and we have defined the 
velocity V = L ■ F, with L = £ -1 denoting a mobility tensor. Here, the Brownian force 
f results from the sum of a large number of collisions of the particle with the surrounding 
fluid, each occurring randomly and independently of the others, so that, according to the 
fluctuation-dissipation theorem, we obtain^: 


(f(0) = o ; (f(f)f(f + T))=2*rc<$(T). 


( 2 ) 
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In general, using non-equilibrium thermodynamics notations, x can be identified with a 
generalized flux, J, while ( • V is a generalized force, F; accordingly, L = is the On- 
sager coefficient relating generalized fluxes and forces through the phenomenological linear 
relation, J = L • F. Therefore, the generalized Langevin equation reads: 


J-L F = J, 


(3) 


where J represents the fluctuating part of the flux, with, 

^J(t)^ = 0; ^J(t) J(t + r)^ = 2kTL5(r). (4) 


Note that, by applying the fluctuation-dissipation theorem, we have tacitly assumed that 
the generalized force is conservative, i.e. F = —For example, when we describe the 
motion of charged particles in electrostatics, J is the electric current (equal to the mean 
particle velocity x per unit charge), while F = — V0o is the electric field induced by an 
electric potential, (j) o; in the same way, J can be a heat flux and F = — VT" 1 the inverse 
temperature gradient. 

Now, following the approach by Feynman and Hibbs^ and applying the central limit 
theorem, we see that the probability of observing a certain Brownian force function f (t) 
consists of the following Gaussian distribution, 


IT [f (£)] oc exp 


[f (£) ■ B (r) • f (t + r)] dtdr 


(5) 


where B is a sort of inverse of the variance of the process, 


B(r) 


(f(f)f(f + r)) 1 


1 

2kT 


C ^(r). 


( 6 ) 


Equations ©, © and © reveals that, as the intensity of the fluctuating force is pro¬ 
portional to the local resistance, particles tend to diffuse from regions of high to regions of 
low fluctuations and, accordingly, a system will tend to move along the paths of minimum 
resistance. 

Now, since f (t) and x (t) are linearly related to one another through the Langevin equa¬ 
tion P), the probability fl [x (£)] that the particle follows the path x(f) is proportional to 
II [f (t)]. Consequently, substituting Eqs. © and © into ©), we obtain: 

i y \ 

~4fcT / ( x -V)-C-(x-V)dt 

\ x (d / 


II [x (t)] = G (x, t |x 0 ) exp 


(7) 








5 


where the normalizing term G is the Jacobian^, depending only on the end points (see the 
end of Chapter UTT1) . 

Finally, the conditional probability II (x, i|xo) that the particle moves from an initial state 
Xo at time t = 0 to a final state x at time t will be equal to the sum of the contributions (j6]) 
of all paths connecting the two events, obtaining 


n(x,f|x 0 ) = G (x, f|x 0 ) 


exp 


£ l x ( r )] 

2 kT 


V [x (t) 


( 8 ) 


where the integral is taken over all paths (note that path integrals can be defined rigorously 
even when the paths x(r) are not continuous functions 21 ') such that x (0) = x 0 and x (t) = x, 
with, 

t 

£ [ x (*)] = J £ [x (r), r] dr , (9) 

o 

and, 

£[x(r),r] = i(x —V) ■ C-(x-V). (10) 

In general, C [x (r), r] is (one half of) the rate of energy dissipation at time r as the 
system follows the path x (r), and £ [x (£)] coincides with (one half of) the energy dissipated 
along that trajectory during the time interval [0 — t]. Due to the obvious analogy with 
analytical mechanics, we could also denote C by Lagrangian time rate and £ by action time 
rate. 

The exponential dependence of the path probability on the action time rate has been 
demonstrated directly by Wang and El Kaabouchiu 22 through the direct simulation of the 
random motion of 10 9 Brownian particles in a linear force field. 

In the particular case of linear velocity fields, the same result was obtained by Onsager 
and Machlup^, starting from the Fokker-Planck equation and following Wiener’s original 
derivation of path integration^. In fact, the rate of energy dissipation (fTOl) coincides, apart 
from the sign, with the Onsager-Machlup function, 


£[x(r),r] = -(s-*j-*,.), (11) 

where, 

S = J ■ F; Tj = |J • ID 1 • J; = |F • L ■ F, (12) 

with J = x, F = £ ■ V and L = denoting the generalized flux, the generalized force and 
the Onsager coefficient, respectively. Here, S is the entropy production rate, while \l >j and 
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T p are Rayleigh’s flux-based and force-based dissipation potentials, depending respectively 
on the rate of change and on the state of the system. 

Considering that 

Sdt = I F • xdt — — f V0 O ■ dx = — A</> = — (</>(x) — </>(x 0 )), (13) 

Jo Jo 

we see that the action time rate (J9]) can also be written in the following form, 

t 

£ [x (t)\ = A (</>) + J £' [x (r), r] dr, (14) 

o 

where 

£' [x (t) , t\ = ix • C • X + • C • v = 'Fj + y F . (15) 

It is worth observing that, for imaginary t = — iu , with u real, these equations just give 
Feynman’s path integral representation of a probability amplitude^. 

III. MINIMUM PATH. 

Among all paths, let us denote by y (r) the one that minimizes £. According to the 
Hamilton-Jacobi formalism of classical mechanics, the momentum time rate p along the 
minimum path can be defined as, 

= C • [y - V (y)]. (16) 

x =y 

Clearly, p has the units of a generalized force and, in fact, it coincides with the Brownian 
force f along the minimum path, i.e., p = f mm , where the subscript ”mw” indicates that 
the quantity is defined along the minimum path. 

Now, define the Hamiltonian time rate H (in fact, H has the units of an energy per unit 
time) as H = p ■ y — £ mm , where £ mm is (one half of) the rate of energy dissipation along 
the minimum path, i.e., 

Umin = £ [y (t) , r] = ^ (y - V) • C • (y - V). (17) 

Thus, we obtain: 

ft = 5 (y + v) • < • (y —v) = R ■ < ■ y - A ■ < ■ v, 




(18) 
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where V = V (y), that is, using Rayleigh’s dissipation potentials defined in Eq. (fT2]l . 

W =(«-.,-<M min • (19) 

Note that, if we minimize the action time rate (TT4|) by defining the momentum time rate 
p' = £ ■ y, at the end we would obtain again the Hamiltonian time rate (|T9j) . 

Here we see that, when t —> —t, then y —> —y, while V(y) —>■ V(y), showing that V. is 
time-invariant along the minimum path—'. Indeed, it is intriguing that the evolution of a 
dissipative system can be described in terms of the conservative motion of a pseudo-system, 
whose ’’energy” "H is constant. 

Expressing in terms of the generalized momenta (fT6jh we have: 

H{ y,p) = P - (V+^C'-p) • (20) 

Now, substituting the expression ()20l) into the the canonical equation, p = —cFH/d\ y, 
the minimum path can be determined explicitly. In particular, in the linear case, when £ is 
constant we obtain: 

p + (VV) • p = 0. (21) 

This is the Euler-Lagrange equation, obtained by minimizing the energy dissipated ((9]) . Note 
that, in the isotropic case, when £ = I, Eq. (12T1) can be rewritten in the following simple 
form, 

y = VI/ + yxB, (22) 

where 

U = ^V 2 ; B = -V x V. (23) 

Similar results were obtained by Wiegler,— who studied the motion of Brownian particles in 
irrotational velocity fields, where B = 0. So, the minimum path describes the trajectory of 
a pseudo-particle of unit mass and unit electric charge immersed in an electric field U and 
a magnetic field B. 

Now, in general, the domain of integration of the path integral is composed of all paths 
whose distance from the minimum path is of order 5 ~ kT/Cy or less, where ( and V are 
typical values of £ and V, respectively. Therefore, expressing any path x (r) as the ’’sum” 
of the minimum path y (r) and a ’’fluctuating” part x (t)2£, 


x (r) = y (r) + x (r), 


(24) 


where x (0) = x (t) = 0, then £ [x (t)] can be expanded formally around y (r) as: 


£ [x (t)] = £ min + -xx: 


d 2 £ 

fix fix 


+ 


(25) 


x =y 


with £ min = £ [y (t ), f], where we have considered that the first order derivative is identically 
zero. Accordingly, we see that, if within distances of O (5) from the minimum path, x can 
be approximated as a linear function, then £ is a quadratic functional, and therefore the 
expansion (|25l) terminates after the second derivative, with the last term being a function of 
x only, and not of y2S. Finally, substituting (l25]i into m-m, with x(r) = y(r), we obtain: 

t 

1 

II (x, t|x 0 ) = G (x, t|x 0 ) exp 


2kT 


^min id~') diT 


where the normalization function G depends only on the end points. 
Subjected to a posteriori verification, we have: 


G (x, t|x 0 ) = exp 


1 

2 


V • (x + V) dt 


Therefore, substituting this expression into Eq. (|26jl . we obtain: 


(26) 


(27) 


II (x, i|x 0 ) = exp 


kT 


— = f (Yj dr 

2 °mzn / ^mm V ; a 1 ? 


(28) 


with, 

dl = \ (y - V) ■ c ■ (y - V) + kTV ■ (y + V), (29) 

and y (t) represents the trajectory (l2l~j) along the minimum path. Note that the last term 
in the RHS of Eq. (|29l) does not enter the minimization process, as its contribution only 
depends on the end points. 

This result shows that under very general conditions the path integral is determined 
exclusively by the minimum path f)2T]l . determining the Boltzmann-like distribution f[28[) . 
When compatible with the end points, the minimum path is obviously y = V. 

Note that here, since the force held is conservative, i.e., F = £ ■ V = — V0o> with 
0o = 0o( x ) denoting a time-independent potential, we obtain: y • £ ■ V = —d<fio/dt. 
Then, the steady state, equilibrium probability distribution of unconstrained systems can 
be obtained considering the reverse minimum path, y = —V, assuming that initially the 












9 


system is at equilibrium, i.e., x (t = 0) = x 0 = 0, with 0 o (O) = 0, while x (t) = x. 
At the end, since £ min = y • V0o = Po, we obtain the usual Boltzmann distribution, 
n(x) = W exp [-0 o (x)/fcT]. 


IV. THE FOKKER-PLANCK EQUATION. 


Now, instead of keeping both the lower and the upper limit of (BHD fixed, let us compare 
the different values of p corresponding to trajectories having the same initial configuration, 
but with variable final configurations, i.e. when Jx(0) = 0 and <5x(£) = Jx. In this case, 
instead of Sp = 0, we obtain Sp = • <5x, and therefore, 

V0 = ^p. (30) 


In the same way, p can be seen as a function of time, i.e. we consider the trajectories 
starting from the same initial configuration and ending at the same final configuration, but 
at a different time. From Eqs. (0 and (BHD, we se that the total derivative of p is: 

dp _ 1 ~(t) 
dt ~ 2 rni 

Note that, as the upper limit of the trajectory is free, the full expression (|29]l of the 
Lagrangian has to be considered. Therefore, expressing the action rate, (p , as a function of 
both spatial and temporal coordinates of the upper limit, i.e., 0 = 0 (x, t), we find: 

dip dp dp . dp 1 

Tt = a + &T x= a + 2 p ' x - 




thus obtaining the following Hamilton-Jacobi equation, 

^ + iw«(x,V0,() = 0. (33) 

Here, U® is the total Hamiltonian time rate, 

^ (f) (x,p,t) = p • x — C® in — p • Qc'-P + v) — 2kT V • Qc^'P + v) , (34) 

with p = 2V0. Therefore, Eq. (|33l) becomes, 

^ + V0 • (C" 1 • V0 + V) = kTW ■ (C 1 • V0 + V) . (35) 
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Now, considering that the conditional probability is related to the generalized potential, 
4>, through II oc exp (— <f>/kT), it is easy to see that the Hamilton-Jacobi equation (135|) is 
identical to the Fokker-Planck equation, 

dn 

_ + V ■ J = 0; J = vn- kT C 1 ■ vn, (36) 

where J is the probability flux. Therefore, we see that the quadratic approximation (1281) . i.e. 
assuming that the minimum path is the only trajectory that contributes the path integral, is 
equivalent to assuming that the Fokker-Planck equation is a valid approximation to describe 
the random process. This result justifies assuming the expression (l27)i for the normalization 
term. 


V. STEADY STATE. 

In this section we consider a system that is kept in a state of non-equilibrium through 
a set of constraints about its generalized forces or fluxes. In this case, considering that the 
time integral in Eq. (j5j) reduces to a simple product, i.e., £ = £t, the rule of least energy 
dissipation becomes: 

£ (x, x) = ^ [x - V (x)] ■ C ■ [x - V (x)] = min (37) 

This leads to Onsager’s principle, 

— £ = (s — T j — j = max, (38) 

where S, j and T p are defined in (1T21) . In addition, considering that at steady state 
d(j)/dt = d(f)/dt = \£mi n i Eq- (1331) becomes: 

£min + ^ W =0, (39) 

so that we obtain: 

p • y = o. (40) 

Consequently, considering that [cfr. Eq. ([16]) ] p = C ■ (y — V), we see that S = 2Tj, so that 
Eq. (1381) can be rewritten as: 


% = (T j — ^ F ) = max, 


(41) 
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where Eq. (jT9|) has been considered. This shows that at steady state a system tends to 
maximize its Hamiltonian time rate. The same result could be obtained from Eq. ([20]) . (I38j) 
and (j39l) . 

Equation (T4TT) expresses a principle of minimum resistance, revealing that if in a process 
the generalized forces, F. are fixed, then the fluxes, J, are maximized, while when the gen¬ 
eralized fluxes are fixed, then the forces are minimized. This confirms that, since from Eqs. 
©, (J5[) and fl6[) the intensity of the fluctuating force is proportional to the local resistance, 
a system will tend to move along the paths of minimum resistance. So, for example, if the 
temperature (or concentration) difference between two regions is fixed, heat (or chemical 
species) will tend to flow across regions of low resistance (or large conductivity), thus max¬ 
imizing the flux. In that case, therefore, the entropy production rate will be maximized 
and, in fact, this behavior satisfies the maximum flux principle. Alternatively, when the 
heat flux is fixed, the rule of minimum resistance leads to minimizing the temperature (or 
concentration) difference, so that entropy production rate will be minimized. 

Studying a similar problem, Martiouchev and Seleznev 2 - showed that the morphologies 
selected during crystallization tend to maximize the entropy production. Similar results 
were obtained by Molin and Mauri^, who simulated the spinodal decomposition of a binary 
mixture confined between two walls that are quenched below the critical temperature. This 
rule was demonstrated in great generality by Dewar^ using nonequilibrium statistical me¬ 
chanics, showing that it applies to all cases where macroscopic fluxes are free to vary under 
imposed constraints. Similar results were obtained by Favretti^. Accordingly, despite its 
name, this law generalizes Onsager’s principle of least energy dissipation- to cases where 
systems are far from equilibrium and, in fact, it agrees with the steepest entropy ascent 
(SEA) model developed by Beretta^. 


VI. UNSTEADY CASES. 

A. Relaxation to equilibrium. 

In this section we consider unconstrained systems that relax towards equilibrium. In 
this case, since at equilibrium p = V = 0, then H = 0. Therefore, considering that 
H is constant along the minimum path, we conclude that during its relaxation towards 
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equilibrium a system satisfies the following relation: 


1-L — 0 —> T j — \Eg, 


(42) 


where Eq. (TT9|) has been considered. 

As a simple example of application, consider a typical Ornstein-Uhlcnbeck process, where 
a Brownian particle, with drag coefficient Qj = (Sij and diffusivity D = kT/(, is immersed in 
a quiescent fluid and is subjected to a linear potential force, (i.e., it is a harmonic oscillator), 
attracting the particle towards the origin, with F = = —CM ■ x, where Mij = Mji. 

Accordingly, we can choose a reference frame where the axes coincide with the principal 
directions of the M matrix, so that = MiSij, with Mi > 0. 

Now, the conditional probability function fl (x, f|0) that describes the motion of this 
Brownian particle is given by Eqs. (I28l) - (l29l) . Here, the minimum path equation (l22jl reduces 
to: 

Vi = Mfyi, (43) 


which, coupled to the conditions r/*(0) = 0 and yi(t) = x i: yields: 

sinli (MjT) 


Vi{r) = Xi 


sinli (Mjt) 


Substituting this result into Eq. (J29]), i.e. C m in — C (y — V) , we obtain: 

M 2 r 2 

Cmin = s i n h 2 '(M t) ^ COsh ^ MiT ) + sinh ( M * r )) 2 + C (r), 


(44) 


(45) 


where C(t) is an irrelevant time-dependent function, independent of the endpoints. Finally, 
from Eq. (12SD we find the following Gaussian distribution^: 


n(x,t|0) 


W (■ t) exp 


-Jp MiX2 i t 1 + COth ( M ^)] > 

i 


(46) 


where W(t) is a normalization factor, which is independent of the endpoints. In particular, 
for long times, t M -1 , this solution tends to the equilibrium distribution, 


H(x) 


W exp 


1 

2 kT 



(47) 


where g t = CM* represents the spring constant along the i-th direction. This is the well- 
known Ornstein-Uhlenbeck solution—, obtained by solving the Fokker-Planck equation; 
clearly, this result is not surprising, since in Eqs. ()35jl and (j36j) we saw that the Fokker- 
Planck equation is equivalent to the Hamilton-Jacobi equation for the mininum path. 









13 


VII. UNIFORM VELOCITY FIELD AND SIMPLE SHEAR FLOW. 


In this section we want to show that the path integral approach can also be applied to 
cases where there is no steady state. The simplest case is the diffusion of a Brownian particle 
in a uniform velocity field, i.e., when V is a constant vector. In this case, the minimum 
path y(r) satisfies the Euler-Lagrange equation f[22lh y = 0, with y(0) = 0 and y (t) = x, 
obtaining: y(r) = xr/t. Thus, the minimum path corresponds to a trajectory with uniform 
velocity, which of course has nothing to do with the ’’real” velocity of a diffusing Brownian 
particle. Now, substituting this result into (I28l) - (l29l) we obtain the well known result: 

(x-Vtf 


n (x, t|0) = W(t) exp 


4 Dt 


(48) 


where W(t) = (AnDt) 3 ^ 2 is a normalization factor, independent of the endpoints. 

A more complex case consists of a Brownian particle immersed in a simple shear flow 
held, Vi = 7 x 2 ; V 2 = V 3 = 0. As shown in Mauri and Haber-22, following the same procedure 
as before, we see that the minimum path equation (l 22 l) reduces to: 


Vi - 72/2 = 0; y 2 + yyi - 7 j/ 2 = 0; y 3 = 0, 


(49) 


which, coupled to the conditions y ?; (0) = 0; y t (t) — x t , yields: 

Vi(t) = Ci (r 3 - 6r) + Car 2 ; y 2 [r) = 3Cir 2 + 2C 2 r; y 3 (r) = Y 3 r/t, (50) 


with t = 77 ", where 


Ci = 


tx 2 — 2x\ 
t (t 2 + 12) ’ 


C 2 = 


(6 — t 2 ) x 2 + Atx 1 


(51) 


and t = 7 1 . Substituting this result into Eqs. 
distribution: 


t (t 2 + 12) 

and (I29jl . we find the following Gaussian 


n (x, t\0) = W ( t ) exp 


37 (2xi - tx 2 ) x\ + 


Xn 


ADt ( t 2 + 12) 
Therefore, the variances of this distribution are: 


ADt 


(x 2 ) = 2 Dt 


1 + 3 (7 a 


; (X 1 X 2 ) — Djt 2 ; (xl) — (x\) — 2 Dt. 


(52) 


(53) 


Here we see that, as expected, the mean freedisplace ment in the how direction grows like 
t 3 . Again, as expected, this result coincides with that by Foister and Van de Ven22 and by 
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Katayama and Terauth^, who obtained it by solving the Fokker-Planck and the Langevin 
equations, respectively. 
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